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1 Introduction 



In the physics of radio pulsars, the stationary structure of radio pulsar magnetospheres has 
been investigated in considerable detail. We mean not the vacuum magnetosphere but a mag- 
netosphere filled with a dense electron-positron plasma. This is because the radio emission 
generated in the magnetosphere by a charged particle flux requires the existence of a dense 
plasma production mechanism for its explanation. This mechanism suggested by Sturrock [l| 
and developed significantly by Ruderman and Sutherland 0] is that an electron-positron pair 
can be efficiently produced by a gamma-ray photon with an energy higher than twice the elec- 
tron rest mass in a strong magnetic field near the surface of a neutron star that is a radio 
pulsar. B « lO''^ G g ff. In turn, euergefc gam,„a-ray photous are emitted by eleetrous 
and positrons as they move in the magnetosphere along magnetic field lines with a significant 
curvature. These photons are called curvature ones. The photon emission and absorption in 
a magnetic field accompanied by the production of electron-positron pairs allow a theory of 
stationary plasma generation in the magnetosphere of a magnetized rotating neutron star to be 
constructed. However, from observations of stationarily operating radio pulsars, it is rather dif- 
ficult to understand what the generation mechanism of radio emission in their magnetospheres 
is and how the plasma is produced. Tracing the dynamics of emission development at various 
frequencies would be very important for understanding the physical processes in the magneto- 
spheres of radio pulsars j5|. In addition, many observations of nonstationarily operating radio 
pulsars have appeared recently. These primarily include the so-called intermittent radio pulsars 
whose radio emission is observed only during a certain time interval that exceeds considerably 
the spin period of the star. Thus, for example, the pulsar PSR B1931+24 "operates" for 5 — 10 
days and then remains silent for 20 — 25 days ^, while the pulsar PSR J1832+0029 "oper- 
ates" for about 300 days and remains silent for about 700 days (see, e.g., the review 0]). The 
spindown, i.e., the energy losses, was measured to be significantly different in the periods of 
operation and silence. Since the radio emission power accounts for a negligible fraction of the 
total rotational energy losses of a neutron star, it would be natural to assume that the silence is 
related to the termination of plasma generation in the magnetosphere. The mechanism of losses 
related to the emission of a magnetodipole wave in a vacuum (silence) can then be distinguished 

from the mechanism of losses related to the emission of a pulsar wind an electron-positron 

plasma flow from the magnetosphere (operation). 

Apart from intermittent pulsars, the group of so-called nulling pulsars has long been 
known, in which no radio emission is observed for some time but not so regularly as in in- 
termittent pulsars and for which the difference in spindown has not yet been measured. The 
nulling fraction exceeds 40% for seven of the 23 pulsars investigated in ^ and reaches 95% for 
PSR J1502-5653 and PSR J1717-4054. 

Yet another group of nonstationary radio sources has been observed recently rotating 

radio transients (RRATs), i.e., sporadically flaring radio sources in which the flare repetition 
phase is retained and for which the corresponding periods typical of normal radio pulsars have 
been measured 9|. There is no doubt that these are also rotating neutron stars. Neverthe- 
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less, the nature of their acti vity is completely unclear, as suggested by the existence of many 
dissimilar models (see, e.g., |lCHl4|). 

In our view, all of the described sources show nonstationary plasma generation processes 
in the magnetospheres of neutron stars. Therefore, it is important to understand how the 
magnetosphere of a rotating magnetized neutron star is filled with plasma, which triggers the 
operation of radio pulsars [5,]. This paper is devoted to an initial investigation of the "ignition" 
of neutron star magnetospheres — the dynamics of magnetosphere filling with electrons and 
positrons produced in the magnetosphere. In the Section 2, we give the electromagnetic fields 
of the inner vacuum magnetosphere and describe the force-free surface where the electric field 
along the magnetic field becomes zero. Sections 3 — 7 are devoted to the dynamics of charged 
particle motion near the force-free surface and their capture onto this surface. In the Section 8, 
we discuss the dynamics of plasma accumulation in the inner magnetosphere and estimate the 
characteristic ignition time. 



2 Structure of the vacuum magnetosphere 



It will be convenient for us to consider the electromagnetic field around the neutron star in 
spherical coordinates r, 9, and ip. We will take the rotation axis of the star specified by the 
direction of its angular velocity vector O as the polar axis. Here, r is the distance from the 
stellar center to a given point, 9 is the polar angle measured from the rotation axis, and (p is the 
azimuthal angle. The electromagnetic field outside the neutron star was found by Deutsch [3] 
to have the following form: the magnetic field 
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(2) 



E^ = k—^ y— - ij cos6'sin6'mSin(v9 - 



Here, 9m and (pm are the polar and azimuthal angles of the magnetic axis specified by the 
direction of the magnetic dipole moment vector m. A; = fi/c is the wave number corresponding 
to the angular frequency of neutron star rotation fi, and R is the stellar radius. The azimuthal 
angle = was chosen in such a way that its value was zero at t = 0. Equations ([T]) and ([2]) 
are valid as long as we consider the electromagnetic field at distances much smaller than the 
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light cylinder radius Rc = c/fl: 



{Qrjcf < 1 



(3) 



where r± is the distance from the rotation axis of the neutron star to the point under consid- 
eration. We are primarily interested in the regions of the inner vacuum magnetosphere where 
an efficient one-photon production of electron-positron pairs is possible. For pulsars with a 
characteristic surface magnetic field of 10^^ G, the characteristic distance at which the pair 
production is still possible is ~ (10—20)/?; therefore, for typical values of the stellar radius, 
i? ~ 10 km, and spin period, P ~ 0.1—1 s, parameter ([3]) is ~ 10~^— 10~^. Under these 
conditions, Eqs. ([1]) and ([2]) yield a fairly accurate result. 

Using Eqs. ([T]) and ([2j) for the electromagnetic field by Deutsch, we can easily derive an 
equation of the force-free surface E ■ B = — the surface at each point of which the longitudinal 
electric field is zero: 



The structure of the force- free surface (jH) can be understood from Fig. 1. It shows the sections 
of the force-free surface by the plane ip — ipm = {0, vr} that passes through the rotation and 
magnetic axes and by the plane <f — = {7r/2,37r/2} that also passes through the rotation 
axis but orthogonally to the preceding plane for various inclinations of the magnetic axis to the 
rotation axis of the neutron star. An idea of the 3D structure of the force-free surface can be 
gotten from Fig. 3 (see below), which shows the particle trajectories lying on this surface. 

For an aligned rotator, the force-free surface is just the equatorial plane 9 = 9' = ti /2. 
For an arbitrary oblique rotator, the force-free surface can be divided into two regions. One 
region consists of two vaulted parts adjacent to the neutron star surface at the points of the 
equator 9 = -n j2 and the magnetic equator 9' = 7r/2. The force- free surface touches the neutron 
star surface at all points of the magnetic equator, which is not observed at the points of the 
ordinary equator. The second region is the surface that starts at the magnetic equator and goes 
to infinity. For an easier association of this surface with the sections by the plane ip — ipm = 
{0, vr}, it is convenient to imagine it as being composed of two open sheets. One edge of each 
sheet is adjacent to the magnetic equator and to the straight line passing through the center of 
the neutron star orthogonally to its rotation and magnetic axes. The other edge of the sheet 
goes to infinity in such a way that the angle 6'" — > 7r/2 as r ^ oo for an arbitrary point lying on 
the force-free surface. At large r, the sheet differs little from the plane 9" = 7t/2. However, at 
relatively small distances from the stellar plane, the difference is noticeable and a dome-shaped 
part ending at the magnetic equator is formed above the plane 9" = tt/2. The two sheets are 
smoothly jointed at the points of the straight line 9 = 9' = 9" = -k /2 mentioned above to form 
a single sheet. 




(4) 



Here, we introduced the angles 9' and 9" as follows: 



cos 9' = cos 9 cos 9m + sin 9 sin 9m cos{ip — ipm), 
cos 9" = — cos 9 sin 9m + sin 9 cos 9m cos{ip — ipm)- 



(5) 
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Below, we will everywhere consider only the values of 6^ lying within the range from 
to 7r/2. This does not limit the generality, because a rotator with an angle x > 7r/2 between the 
vectors m and fl is equivalent to a rotator with 6'm = 7r — x<7r/2; —Q should be substituted 
for Q in all formulas. 



3 The motion of charged particles 

To investigate the motion of particles in a vacuum magnetosphere, we will use the classical 
Dirac-Lorentz equation 

o r 1 

(6) 
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where = (ct, r ) is a contravariant 4- vector containing the time t and coordinates r = 
{x,y,z)^ of the particle in the laboratory frame of reference (T is the transposition symbol), 
rrie is the particle mass, e is the particle charge, and c is the speed of light. The dot over the 
4- vector denotes differentiation with respect to the proper time r of the particle, i.e., the time 
in the comoving frame of reference. The differentials of the time dt in the laboratory frame of 
reference and the proper time dr of the particle are related by the relation dr = dt/'-f, where 7 
is the particle Lorentz factor. The 4-force acting to the particle is defined as = eF^^Xk/c, 
where 

F,k = dAk/dx' - dAi/dx'' 
is the electromagnetic field tensor and the definition of the 4-potential 

A' = {A^,A\A\AY = {<P,^^Y 



includes the standard scalar, 0, and vector. A, electromagnetic potentials. 

Let us pass to dimensionless variables. We will measure the strength of the electric 
and magnetic fields in units of the so-called critical field Be = mlc^/eh ~ 4.4 x 10^^ G, the 
particle velocity in units of the speed of light c, the particle charge in units of the positron 
charge e, the particle mass in units of the electron mass rUe, the particle energy in units of 
the electron rest mass meC^, all distances in units of the electron Compton wavelength X = 
h/rrieC ~ 3.86 x 10^^^ cm, and all times in units of X/c. We will immediately note that in these 
units, 1000 km ^ 2.6 x 10^^ and 1 s ^ 7.8 x 10^°. 

After the reduction to dimensionless form described above and the separation of the 
scalar and vector components of the 4- vector x\ the Dirac-Lorentz equation ([6]) is reduced to 
the system of equations 



d'j 
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± 



E- v(v-E) + v X B 



(7) 
(8) 



where v is the particle velocity and a = e^/hc ~ 1/137 is the fine-structure constant. When 
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the motions of positrons and electrons are investigated, the plus and minus signs, respectively, 
should be taken in the system of equations ([7j) and ([H])- For convenience, we will write the 
plus sign everywhere below, separately stipulating in important cases what will change if not 
positrons but electrons are considered. 

A charged particle produced in the magnetosphere, be it an electron or a positron, will 
be subjected to an electric field and will accelerate. It is easy to estimate the time in which 
the particle motion will become relativistic, r^e/ ~ From the general form of Eqs. ([1]) 

and ([2]), we see how the electric field strength E is related to the magnetic field strength B, 
with the relation Eg^rf ~ RBsurf/Rc holding for the fields on the stellar surface. Since the 
surface electric field for a typical surface magnetic field Bgurf ~ 0.01—0.10 and a ratio of the 
neutron star radius to the light cylinder radius R/Rc ~ 10^"^— 10^'^ is Esurf ~ 10^^—10"^, the 
particle transition time to the relativistic regime is r^ei ~ 10^—10^. This means that the particle 
reaches a relativistic velocity in a time of lO^^'' — 10^^^ s and subsequently its motion may be 
considered as ultrarelativistic. Obviously, in this case, the particle will travel a distance of no 
more than 10^—10^ Compton wavelengths. 

Subsequently, the particle will continue to accelerate already in the ultrarelativistic 
regime. As is well known, in general, the particle will move along curved magnetic field lines. 
Acquiring an increasingly high energy, the particle will begin to emit the so-called curvature 
photons. Intense curvature radiation leads to the loss of particle energy. The emerging force of 
radiative friction is taken into account by the Dirac-Lorentz equation. 

The intensity of curvature radiation increases with Lorentz factor 7. Therefore, in the 
long run, having reached some maximum Lorentz factor 70, the particle subsequently will not 
undergo any acceleration, since the entire energy gained by it through the action of the electric 
field will be lost through the emission of curvature photons. To determine 70, we should find a 
stationary solution of Eq. ([?]) by setting d'y/dt = d'^'j/dt'^ = 0: 



Equation ([9]) shows that in a stationary state, the entire work of the electric field on the particle 
completely transforms into the energy of curvature radiation. In the case under consideration, 
the particle velocity vector has a unit length, i.e., v = 1. Therefore, d\/dt = n/p, where n is 
the vector of the principal normal to the particle trajectory and p is the radius of curvature of 
the trajectory. The maximum particle Lorentz factor takes the form 



For a characteristic longitudinal electric field E\\ ~ 10~^ and a radius of curvature of the particle 
trajectory p R 2.6 x 10^^, the maximum particle Lorentz factor is 70 ~ 6 x 10"^. 

Let us estimate the time it takes for the particle to reach the maximum Lorentz factor 
7o — the time of its full acceleration and its transition to the quasi-stationary regime of motion 




(9) 




(10) 
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defined by the balance condition ([9]). The sought-for time is r^j ~ 7o/-E'||. Thus, the particle 
reaches its maximum energy in a time r^t ~ 10^^ (10~^ s in dimensional units) and subsequently 
moves in such a way that the work of the electric field on the particle per unit time is equal 
to the total intensity of curvature radiation. In this case, the particle travels a distance of 
~ 10^^ (of the orders of several tens of centimeters in dimensional units), which is much smaller 
than R. Therefore, the assumption that the electric field does not change in the time r^j of full 
particle acceleration is valid. Hence, the electrons and positrons may be assumed to accelerate 
almost instantaneously at the point where the electron-positron pair was produced. 

However, since both the electric field strength and the radius of curvature of the trajec- 
tory change as the particle moves, the Lorentz factor 70 also changes with time. This means 
that the particle will adjust itself, acquiring energy due to the work of the electric field on it if 
7o increases along the trajectory and losing energy through the emission of curvature photons 
if 7o decreases along the particle trajectory. Now, it is important to determine the rate of this 
adjustment of the particle energy. If this rate exceeds considerably the rate of change in 70 
along the trajectory, then the particle Lorentz factor may be assumed to be determined by 
the coordinates of the point at which the particle is located, because the value of 70 itself at 
some fixed instant of time depends only on the coordinates of the point under consideration 
but under no circumstances on the particle velocity. 

Let us find the law of change in the particle Lorentz factor 7 as its stationary value of 70 is 
approached. For this purpose, we should use Eq. ([7]) by representing the particle Lorentz factor 
as the sum of its stationary value, 70, and some deviation from it, ^7. The deviation ^7 should 
be considered small, i.e., ^7 <^ 70. Linearizing Eq. ([7]) and eliminating the self-accelerating 
solution, we obtain ^7 = 57^ exp(— t/ro). Thus, when the particle Lorentz factor deviates from 
its stationary value, 70, by 57i, it will approach 70 exponentially with a decay time constant 
To = 3p^/8a7o. For a characteristic radius of curvature of the trajectory p ~ 10^^ and Lorentz 
factor 7o ~ 10^—10^, the decay time is tq ~ 10^^ — 10^^ (10^^° — 10^^ s in dimensional units). In 
time To, the particle will travel a distance Iq = tq of ~ 10^^ — 10^^ Compton wavelengths (i.e., 
from centimeters to tens of meters). Now, recall that the characteristic distances at which the 
electric field and the radius of curvature of the magnetic field lines and, hence, the stationary 
particle Lorentz factor 70 change are of the order of R. Since /q ^ -R, the particles will adjust 
themselves to a change in 70 as they move in the magnetosphere. Thus, the Lorentz factor of 
the particle may be assumed to be completely determined by its coordinates. 

So far, we have discussed the particle energetics using only Eq. ([7]). Let us now in- 
vestigate the second equation ([8j) in more detail and examine how this equation defines the 
particle trajectory. Let us estimate the terms that arise when the self-action of charged parti- 
cles is taken into account by assuming that the time r^f has elapsed after the generation and 
the particle has already fully accelerated and passed to the quasi-stationary regime of motion 
defined by condition ([9]). The first term in Eq. ([8]) containing the square brackets is equal in 
order of magnitude to a'-^^/R^ . It is small compared to the second term, because the condition 
u'-Jq/R^ <^ E is met. This can be verified by taking typical values of the particle Lorentz 
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factor, 7o ~ 10^, the stellar radius, R ~ 10^"^, and the electric field, E ~ 10"''. However, the 
satisfaction of this condition is not a sufficient reason to discard the self-action terms. Indeed, 
the term on the left-hand side of Eq. ^ is equal in order of magnitude to 70/-R ~ 10"^ and is 
also much smaller than E. For a justified discarding of the terms that allow for the effect of 
the field generated by charged particles on their proper motion, it is necessary that these terms 
be small compared to 70/-R, ct'-yo/R ^ 1. We see that this condition is always met; therefore, 
the motion of particles in the vacuum magnetosphere of a neutron star will be described by the 
equation 

dv 

7— = E- v(v-E) +v X B. (11) 

Thus, Eq. ( TTT]) describes the proper motion of the particle, while its energetics is defined by 
Eq. ([7]), which is reduced to equality ffTOl) in the quasi- stationary case.. 

The total particle velocity v will be rewritten as v = b + Vg + v^, where b = B/i? is a 
unit tangent vector to the magnetic field line, Ve = E x B/i?^ is the electric drift velocity, and 
Vc = 7b X dh/dt/B is the centrifugal drift velocity. 

Let us estimate the drift velocities Vg and v^. in order of magnitude. The electric 
drift velocity is equal to the ratio of the electric and magnetic field strengths, Ve ~ E/B, and 
is R/Rc ~ 10-^ in order of magnitude. For 7 ~ 10®, p ~ 10^^ and B ~ 0.01-0.10, the 
centrifugal drift velocity is Vc ~ 10~®— 10"''. Thus, the centrifugal drift velocity is a quantity of 
the next order of smallness compared to the electric drift velocity; therefore, we will disregard 
the centrifugal drift. So, we will assume that v = b + Vg to terms o{kr) — quantities of the next 
order of smallness (quadratic or higher order) in kr ~ R/Rc. 

The electromagnetic field described by Eqs. ([T]) and ([2]) depends periodically on time. 
Let t, r, 6, and if be the spherical coordinates in the laboratory frame of reference. Let us pass 
to new coordinates t', r', and (f' using the formulas t' = t, r' = r, -d' = 6, and ip' = (f — Qt, 
where t', r', d' , (p' are the time and the coordinates in a rotating frame of reference. In these 
coordinates, the electromagnetic field depends only on r', d\ and if' but not on t' . Below, we 
will omit the primes on the variables in the cases where it is clear that the analysis is performed 
in the rotating frame of reference. The velocity and acceleration transformations when passing 
to the rotating frame of reference have a standard form: 

V = V + \tr, 

— = — + 2nxv' + rjx vt„ 

dt dt' 

where v' is the relative velocity and \tr = x r is the translational velocity. Formulas (fT2|) are 
valid ii VL <^Vt^. This condition is equivalent to the requirement that P <^ 1, which is always 
met, since P ~ 10~^^ in order of magnitude. Substituting Eqs. (1121) into Eq. flTTl) yields the 
following equation of particle motion in the rotating frame of reference: 

7— = E + v,, xB-v'(v'-E) + v'xB. (13) 
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Here, E and B are the electric and magnetic fields defined by Eqs. ([2]) and ([T]), respectively, 
but now, after the change of variables described above, they depend on the coordinates r', d' , 
and if'. In deriving Eq. fll3p . we neglected the Coriolis, 20 x v', and translational, f2 x v^j., 
accelerations compared to the relative acceleration dV /dt' . This is possible, because the relative 
acceleration is mainly the axipetal acceleration that arises from the motion of particles along 
curved magnetic field lines. Therefore, \dv'/dt'\ ~ l/i?andis~ 10^^^. The Coriolis acceleration 
is |2ri X v'l ~ 1/-Rc in order of magnitude, which is ~ 10"^^. The translational acceleration 
is even lower, because the translational velocity vtr is definitely lower than f' ~ 1 in the 
magnetospheric regions under consideration, r < lOi? ^ Rc-, where an efficient one-photon 
production of pairs is possible. The term 7 x (2v' + v^^) being discarded, which would be 
present on the left-hand side of Eq. (fT3i) . is also small compared to the electric field strength, 
E ~ 10~^, since it is 7/-Rc ~ 10^^^ for 7 ~ 10*. The Lorentz factor 7 in Eq. (I13p refers to the 
laboratory frame of reference. 

Below, we will investigate the motion of particles in the rotating frame of reference. We 
are interested in the question of whether the regions in which the accumulation of a primary 
plasma is possible exist in the magnetosphere. Let us find the equilibrium positions — the 
points in moving to any of which a charged particle will remain at it for an unlimited time. 
The coordinates of the equilibrium positions are determined by the conditions for the particle 
velocity and acceleration being zero in the rotating frame of reference: v' = dw' /dt' = 0. 
Substituting this condition into Eq. (fT3|) yields the equation = 0, where the effective 

electric field is E^-^-'^ = E + v^^ x B. Its components, which can be easily obtained from Eqs. ([1]) 
and ([2]) using the expression for v^^, in spherical coordinates are 



where Br, Bg, and B^, are defined by Eq. The condition E"^^ = gives the following set 
of equilibrium points: 




(14) 




(equator) , 



(magnetic equator). 



(15) 



r+ _ 3 + cos 9., 
^ ~ 1 -cos^, 



m 



m 



(open sheets). 



(16) 





3 — cos 6, 



1 + cos 6. 



m 



m 



(vaults) 



(17) 
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We see that all equilibrium points are on the force-free surface (jl]) and are separated into two 
groups: nonisolated and isolated. All points at the equator and the magnetic equator are 
nonisolated equilibrium points (see (ITB]) ). Four more isolated equilibrium points are added to 
these two circumferences: two points with the coordinates given by Eqs. f|T6l) are located on the 
open sheets of the force-free surface and the other two points defined by Eqs. (1171) are located 
in its vaulted parts whose boundary is adjacent to the equator and the magnetic equator. All 
four points lie in the plane passing through the rotation and magnetic axes. 

The question arises as to whether the equilibrium positions found are stable. In addition, 
the behavior of the particle trajectories near the force-free surface in general and near the 
equilibrium positions in particular is not yet clear. We will give an answer to this question 
after a detailed study of the particle capture by the force-free surface. 



4 Oscillations of charged particles near the force-free 
surface 

We will quantitatively investigate the motion of charged particles near the force-free surface in 
the rotating frame of reference using Eq. (I13p . Let us choose some point fq on the force- free 
surface and decompose the electric field E, the magnetic field B, and the effective electric field 
Ei^f^ introduced above near this point as 

E = Eo + (x' ■ V)Eo, 

B = Bo + (x' • V)Bo, (18) 
E^^^ = E;;^^ + (x' ■ V)El^^, 

where Eq, Bq, and Eq-^-'^ are the field strengths at point fq; and x' = r — fq is the distance 
from point r at which we are interested in the field strengths to point tq. Decompositions (ITSj) 
are valid, because we consider the particle motion in a small neighborhood of point — at 
distances much smaller than the characteristic distances of change in the fields, so that x' <^ R. 
Before seeking for the complete solution of Eq. (fT3l) . let us first find its particular solution. We 
will require that this particular solution describes the particle motion with a constant velocity, 
Vq = const. For this solution to exist, two conditions must be met simultaneously: 

(v[,-V)E^^^ + v^x (v[,- V)Bo = 0. 

The first equation of system f|T9l) uniquely defines the velocity component orthogonal 
to the magnetic field direction, = Eq-^-^ x Bq/Bq, with the longitudinal component v'n still 
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being arbitrary: Vq = fyb + v^. The components of the velocity are 

1 n 



v'^Q = Qr 



cos 9br + (1 J ] sin 9 6, 



?2s 



(20) 

(^1 - — j sin 6 (^-6^ + 62 j + ^os ^ 



where br, be, and 6;^ are the components of the unit vector b. 

The longitudinal component v'^^ of the velocity Vq is uniquely defined by the second 
equation of system f[T^ . It is easy to obtain the relation 

v[, ■ V(Eo ■ Bo) = 0. (21) 

Since the gradient V(Eo ■ Bq) is directed along the normal to the force- free surface 
E ■ B = 0, the velocity Vq lies in the tangent plane drawn through point tq of the force-free 
surface. 

Thus, the particular solution that describes the particle motion with a constant velocity 
exists. The transverse velocity is given by Eqs. fl2U]) . while the longitudinal velocity can be 
obtained from condition fl2Tl) and is 



vl ■ V(Eo ■ Bp) 
■ b ■ V(Eo ■ Bo) 



We see that the particle velocity Vq lies in the tangent plane. Therefore, the particle 
cannot leave the force-free surface and its trajectory lies entirely on this surface. The velocity 
will change from point to point and is defined by the same expression Vq = f yb + v^, but the 
field strengths at the point where the particle is located at a given instant of time should be 
taken every time as Eq-^-'^ and Bo. 

The complete solution of Eq. (|T3l) can be represented as 

x' = x[, + x;, V = V, + V„ (22) 

where Xq and Vq are the adiabatic solution found above. No special constraints, except the 
smallness of x'^ compared to the characteristic scales of change in the electromagnetic field 
equal to R, are imposed on x'^^ and v'^ a priori. In particular, the velocity v'^^ can be close to 
the speed of light. Substituting equalities fl22|) into f|T3|) yields the equation 



7 



^ = (x; . v)E^^^ + X (x; ■ v)Bo + v; X B - v^v; ■ e), (23) 



where E and B are defined by Eqs. fllSp . while the coordinate x' in them is defined by f l22p . In 
deriving Eq. fl23l) . we used the condition v'q <^ 1. The term v'^(v']^ ■ E) is retained, because, in 
general, the velocity v'^ can be of the order of unity; therefore, the term under consideration 
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can be equal to E in order of magnitude. 

To derive an equation for describing the oscillations of the particle with no special 
constraints imposed on its velocity, we must use Eq. (I23p by assuming that v'^ || B and write 
the derived equation in projection onto the magnetic field direction: 

7^ = -^M - (24) 

Here, uj is the frequency of nonrelativistic oscillations defined by the formula 

2 b ■ V(Eo ■ Bo) 

u = . (25) 

Bo 

Let us first consider the nonrelativistic case of particle oscillations. Since the term 
{v[yE\\ at f ^ ^ 1 is negligibly small and 7^1, Eq. flM|) transforms to an ordinary equation 
of nonrelativistic oscillations with frequency u. Let us find the applicability criterion for the 
nonrelativistic approximation. As we see from Eq. (1251) . the characteristic oscillation frequency 
is cij ~ \/VtB or, equivalently, uj ~ B /Rc. For a characteristic light cylinder radius Rc ~ 10^° 
(for pulsars with a period P ~ 1 s) and magnetic fields B ~ 0.01—0.10, the frequency of nonrel- 
ativistic oscillations is ~ 10^^^ — 10~^° (corresponding to a frequency v = uj/2tt ~ 1—10 GHz 
in dimensional units). Over the oscillation period, the particle will traverse a distance that 
is definitely smaller than 1/z/ and the maximum amplitude of nonrelativistic oscillations is 
Inro ~ l/ct;. The oscillatory motion will be nonrelativistic if x'^ <^ Inro- For the frequencies uj 
found above, we have Inro ~ 10^*^—10^^ (0.1 — 1.0 cm in dimensional units). We see that the 
oscillations of a charged particle near the force-free surface with an amplitude exceeding 1 cm 
are definitely relativistic. This suggests that considering the relativistic case of oscillations that 
is actually realized is of greatest interest. 

Let now the velocity v'^ be not low. We must then use the condition x'^ <^ R and, 
according to Eq. fl25|) . write the longitudinal electric field E\\ = —lJ^x\. Taking into account 



the equality 7 ~ l/yl ~ ("^i)^' obtain an equation for the relativistic oscillations of a 
charged particle near the force- free surface: 

-r- = ^A- (26) 

dt' 73 1 ^ ' 



This equation has the first integral 



C^7 + ^. (27) 



This first integral is equal to the maximum Lorentz factor 7max that a charged particle reaches 
when passing through the force-free surface, i.e., C = 7max5 when x[ = 0. The first integral C 
can also be expressed in terms of the oscillation amplitude A: C = 1 + 

For the convenience of our subsequent analysis, will make the following changes of vari- 
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ables: 



sm 00 = -, sm = -j= — 



K 



C-l C-1 



c + r V c 



(28) 



If we introduce a function 



i?(0, k) = E{(f), k) — COS^0oF(0, k) 

then the relation between and t is given by the equation 



t = ^^(M. (29) 

uj sin 00 

Here, -F(0, k) and -E(0, k) are the elhptic integrals of the first and second kinds, respectively. 

The oscillation period can be accurately determined at once by noting that = 
corresponds to an equilibrium position, while = 7r/2 corresponds to a turning point. Since 
the particle travels the distance from the equilibrium position x[ = to the turning point 
x[ = A in a. quarter of the period. 



UJ sin 00 ' 



where R(k) = E(/t) — cos^0oK(fi;), while K(k) = F(7r/2, k) and E(k) = E{'it/2,k) are the 
complete elliptic integrals of the first and second kinds. 

The exact solution of the oscillation equation (12^ can be written as 

x[ = AsmQ( y (30) 



v/2(C+l) 

Here, we introduced a function Q{z) that is the inverse of the function R{(j),K,), so that z = 
R{Q{z), k) for any real number 2; G M. This can be done, because the function -R(0, k) strictly 
increases and is continuously differentiable with respect to the variable on the entire real M 
axis at ^ K < 1, with i?(M, k) = M. When these conditions are met, the inverse function Q{z) 
exists and is also a single-valued, strictly increasing and continuously differentiable function on 
the M axis, so that Q(M) = M. Its derivative dQ/dz nowhere becomes neither zero nor infinite, 
since this is true for the partial derivative dR{(f), K,)/d(j) at all C from the interval 1 ^ C < 00 of 
interest to us (this corresponds to any physically possible values of the particle Lorentz factor 

Tmax ) • 

Equation (130|) exhausts the problem of charged particle oscillations near the force-free 
surface, provided that there are no energy losses through radiation. The asymptotic limits of 
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the function Q{z) at values of k close to zero and unity are 



2z, K = 0, 

arcsin(z — 2 h(z)) + TT 11(2;), k ^ 1, 



(31) 



where the function h{z) taking on integer values is defined by the formula 

h(z) 



Here, the integer part of a real number y is denoted by [y\. No constraints are imposed on z 
in Eq. fl31l) . In the nonrelativistic case of oscillations, we must use the asymptotic limit of the 
function Q{z) at /t = 0. As we see from Eq. fl30|) . the equality z = ut/2 holds at C ~ 1. We 
then immediately obtain the equation of harmonic oscillations x'l = Asm{ujt) with frequency 
u. In the ultrarelativistic case, C ^ 1, which corresponds to k — 1, and substituting the 
corresponding asymptotic limit of Q{z) fl^ into Eq. (150]) gives 



Ai-1 



2Hz)) 



(32) 



where z = t/A. This equation describes saw-tooth oscillations with amplitude A and period 
T = AA. 



5 Energy losses during oscillations 

During its oscillations near the force-free surface, a charged particle will lose its energy through 
radiation. Let us investigate the time dependence of the rate of decrease in particle energy. 
The particle energy is characterized by the first integral C found above (see fl271) ). The change 
in energy in the ultrarelativistic case where C ^ 1 is of greatest interest to us, since in the 
nonrelativistic case, C ~ 1 and is almost constant. 

Let us turn to the Dirac-Lorentz equation ([7]). Using the fact that ^ 1, we obtain 



d£_2 
~dt ~ 3"^ 



(33) 



The rate of decrease in C will not be constant on timescales of the order of the period of particle 
motion, because it depends on the specific particle Lorentz factor 7, whose value, obviously, 
changes significantly in a quarter of the oscillation period from 1 at the turning point x'^ = Aio 
C ^ 1 at x'^ = 0. Nevertheless, this will not be an obstacle to solving the problem of change in 
particle energy if the value of C itself changes little over the oscillation period, dC/dt <^ C/T. 
When this adiabatic approximation holds, Eq. (l33i) can be carefully averaged over the oscillation 
period. In this case, we must average the expressions dependent on the particle coordinate x[ 
not formally but using the explicit time dependence fl52]) of the coordinate of the oscillating 
particle x[{t), which is an exact one at C ^ L 
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After its averaging, the differential equation for C transforms to 

dC 4 2^ r 64 
dt 9 



105 cj2p^ 



(34) 



In deriving Eq. fl34p .we neglected the second derivative of the energy C with respect to time. 
This action can be justified rigorously only a posteriori, because this requires the knowledge 
of some facts about the structure of the solution C{t). Since this justification, of little interest 
compared to the main results, is cumbersome, we do not provide it here. With the above 
stipulations, the exact solution of Eq. (IMj) is 

e{t)=[{l+eo')e''/^^-l]-'^\ (35) 

where e{t) = C{t)/ Ccurv is the normalized first integral, with Eq = e{0) being its value at the 
initial time t = 0. Here, we also introduced the quantity 

1 nc; 1/3 

Ccurv = ^^M'/' ^ 1.2 M'/' (36) 

and the decay time constant = 9/4ac<j^. The physical meaning of the constants Ccurv and 
can be easily understood by investigating the behavior of solution fl35p at various initial values 
of Eq. Let first Eq <^1. The solution of Eq. flMl) is then given by the formula 



Eit) = Eoe-'/^' (Co <^ Ccurv). (37) 



We see that the case of Cq ^ Ccurv corresponds to discarding the term in brackets in Eq. ([3i 
proportional to {C/CcurvY- However, this means that the energy losses by a charged particle 
through the emission of curvature photons by it are small compared to the losses due to the 
presence of a term proportional to 'jd'^'j/dt'^ in Eq. (155]) . These energy losses can be associ- 
ated with the bremsstrahlung losses attributable to a pure oscillatory particle motion. They 
are nonzero only if the particle energy changes with time. In contrast to them, the curvature 
losses depend only on the radius of curvature of the particle trajectory and the particle en- 
ergy. Obviously, the intensity of the curvature losses does not depend on the time derivatives 
of the charged particle energy. Thus, at Cq <^ Ccurv the particle energy decreases exponen- 
tially with the decay time constant r^, while the losses themselves are attributable mainly to 
bremsstrahlung. 

Let now Eq ^ 1. Obviously, in this case, the behavior of C on relatively short time 
scales t is of greatest interest, because C decreases with time and becomes smaller than Ccurv 
on fairly long time scales, so that the decay becomes exponential. On short time scales t, we 
have 

e{t) =Eo(l- — ) , Tp = ^, 

\ ^pj 4 (38) 

Co ^ Ccurvj t <^ Tp. 
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We see that if Co ^ Ccurv, the particle energy decreases hnearly on time scales much shorter 
than Tp. However, when t becomes comparable to Tp, the decay becomes a power-law one: 

Here, the meaning of the time Tp is clarified as the time in which the initially linear decrease in 
charged particle energy becomes a power-law one. When a time of the order of is reached, 
the dependence changes again and becomes exponential: 

£(t) = e-*/^^ Co>C,,™, t>Td. (39) 

In contrast to the previously considered case of (137|) . at Co ^ Ccurv, the energy losses 
of the particle are attributable mainly to the emission of curvature photons by it. This pro- 
cess dominates on time scales t < t^ and has a power-law pattern. At t > r^^, the particle 
energy decay becomes exponential, indicative of the particle transition to a regime where the 
bremsstrahlung related to the acceleration that the particle undergoes due to its oscillatory 
motion makes a major contribution to the energy losses. Clearly, the transition itself occurs on 
time scales of the order of t^. However, we will make a certain refinement by finding the time 

Tcurv at which C = Ccurv, 1-6., e{Tcurv) = 1, USiug Eq. ([35]), Tcurv = In 2x^/3 ^ 0.23 Trf. 

It should be noted that Tcurv is completely independent of the initial particle energy Eq. 
This suggests that a charged particle with an initial energy Eq ^ 1, even if very high, will 
lose the bulk of it in time Tcurv Subsequently, the particle will have an energy £ = 1 at time 
t = Tcurv and an ordinary exponential decay will then take place with the time constant t^. It 
is not surprising that the dependence (see ( !39l) ) on time scales t > t^ does not contain the 
initial particle energy either. This is because Tcurv is approximately a factor of 4 shorter than 
Td, while E becomes equal to unity in time Tcurv, irrespective of Eq. On long time scales, the 
relation E{Tcurv) = 1 begins to act as the initial condition for the subsequent exponential decay, 
as is shown by Eq. fl39|) . 

All of the results that we obtained by averaging Eq. (l33ll are valid when the condition 
of the adiabatic approximation is met. As follows from Eq. this condition gives an upper 
limit on the charged particle oscillation amplitude: 

A « A_ = (^) uj-'I'p^l' ^ 2.6 uj-'l''p'l\ (40) 
y iZoG; J 

The condition of the adiabatic approximation can also be interpreted as follows. The oscillation 
period must be definitely shorter than the characteristic charged particle energy decay time. 
As such, we must take the fastest time of change in energy that, in our case, is the transition 
time Tp to a power-law decay. Indeed, if we write the condition T Tp, then the upper limit 
on the oscillation amplitude following from it will closely coincide with condition (HUl) . 

Finally, let us consider the energy losses of a charged particle in the case of nonrelativistic 
oscillations that is realized when the oscillation amplitude eventually becomes smaller than Inro- 
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In this case, the time dependence of the coordinate is 

x\ = Aoe"*/^""" cosut, 

where Aq is the initial amphtude and the decay time constant is Tnro = S/au^. As would be 
expected, it is equal in order of magnitude to the decay time constant for the ultrarelativistic 
case of oscillations. 

Let us now estimate all of the quantities introduced above that characterize the time 
dependence of the energy of a charged particle oscillating near the force-free surface in order of 
magnitude. First, let us estimate the energy Ccurv, Ccutv ~ {BB? / R^^^^ . For a characteristic 
surface magnetic field B ~ 0.01 — 0.10, a neutron star radius R ~ 10^'', and a light cylinder 
radius Rc ~ 10^^—10^° (we will use these values everywhere for our estimations below), Ccurv 
is ~ 10^. We can also introduce the corresponding amplitude 

Acurv = -\/'2Ccurv/^ ~ {RRc/ B)^^^ 

for the same magnetic field, neutron star radius, and light cylinder radius: Acurv ~ 10^^ — 10^^ 
(approximately ~ 1 m in dimensional units). 

The characteristic exponential decay time constant, ~ Rc/aB, is 10^^ — 10^^ (~ 
10—10'^ s in dimensional units). The time Tcurv is also of the same order of magnitude. Note 
that the time exceeds the characteristic pulsar spin periods P. Nevertheless, if the fields 
become comparable to critical fields i? ~ 1, which is true for magnetars, then the time can be 
even shorter than the neutron star spin period. Note, incidentally, that the ratio r^/P depends 
only on the surface magnetic field B but does not depend on P. 

The order of magnitude of Tp is not fixed so clearly, because it depends on the initial 
particle energy ^o- However, the range in which Tp lies can be determined easily. The upper 
boundary is given by the condition Eq ^ 1 or, equivalently, Co ^ Ccurv, when it makes sense to 
introduce the concept of time Tp. This means that, in any case, Tp <^ t^. We will determine the 
lower boundary by taking the maximum possible energy Cq. Obviously, this value is definitely 
lower than the maximum achievable Lorentz factor 70 that a particle far from the force-free 
surface will have. Using the fact that 70 ~ 10^ and assuming that Eq ~ ^qI Ccurv ~ 10^, we will 
find that the inequality Tp ^ lO^^^r^ necessarily holds. This lower limit was estimated with a 
margin, since, in reality, an oscillating particle cannot have an energy Ccurv ~ 10*. The reason 
is that even if the particle has such an energy far from the force-free surface, its energy after 
its capture will be slightly lower, since the particle loses part of its energy as it approaches the 
force-free surface. 

Let us also find the order of magnitude of the upper limit Amax on the oscillation 
amplitude A when the adiabatic approximation is still applicable. Using Eq. fj40l) for this 

1 II 

purpose, we will obtain A 

max 

{leRllB^Y'. Thus, A 

max ~ 10^"^ — 10^^ (of the order of several 
tens of meters in dimensional units). An oscillation period T^ax ~ A corresponds to this 
amplitude, which corresponds to Tmax ~ 10^*^— 10^'' s and a frequency v ~ 10—100 MHz in 



17 



dimensional units. Using the amplitude Amax, we easily find the maximum oscillation energy 

1 /7 

C'max ~ (BR'^/Rc) , which is Cmax ~ (3 X 10^)— 10^. This allows us to refine the lower limit for 
the time Tp. For this purpose, we should take a ratio Craa^/Ccurv ~ 10^ — 10^ as the maximum 
energy Eq. We will then obtain a more realistic lower limit, Tp ^ lO^^Trf. 

6 The capture of charged particles 

Above, we found that a charged particle (an electron or a positron) far from the force-free 
surface moves in such a way that the Lorentz factor of the particle is determined only by 
its coordinates and is 70 (see ffTOl) ). However, near the force-free surface, the electric field is 
weak and the particle energy adjustment time tq is comparable to the characteristic distance 
(in dimensionless units) of change in field — in the case under consideration, the distance to 
the force-free surface. This distance Ic can be determined self-consistently from the equation 
Tq = Ic] we should express the longitudinal electric field E\\ on which tq depends in terms of the 
same distance Ic measured from the force-free surface along the magnetic field line along which 
the particle moves. We will use the linear approximation by assuming that Ic is small enough: 
E\\{1) = —u'^l. Here, we formally introduced the quantity u'^ = —dE\\/dl; the derivative dE\\/dl 
should be taken at the point of intersection of the magnetic field line along which the particle 
moves with the force-free surface. It is easy to verify that uj is the frequency of nonrelativistic 
particle oscillations and is defined by Eq. (125|) . Subsequently, we immediately obtain the capture 
length Ic- 



Note that the capture length coincides, to a coefficient, with the maximum oscillation 
amplitude A^ax (HO]) defined by the condition of the adiabatic approximation. This gives an 
insight into how the transition from the ultrarelativistic quasi-stationary motion of a particle 
far from the force-free surface defined by the energy balance condition ([9]) to the oscillatory 
motion of the particle near the force-free surface occurs. For this purpose, we can use qualitative 
considerations. The capture length Ic is almost a factor of 3 smaller than the amplitude v4max- 
This gives us reason to believe that when captured, the particle passes almost immediately to 
the regime of adiabatic oscillations. To verify this fact, let us estimate the capture amplitude 
Ac — the maximum distance to which the particle will be defiected from the force-free surface 
when passing through it. It is implied that the charged particle was initially produced far 
from the force-free surface, started moving toward it, and crosses this surface for the first time. 
Consider the particle that has not yet crossed the force- free surface at distance Ic from it. On the 
one hand, we can estimate the particle Lorentz factor 7^ at a given point from Eq. (ITOl) . using the 
relation E\\{1) = —u'^l, and setting / = Ic, 7c = 4^^^/^. On the other hand, we can calculate the 
first integral flTH) by setting 7 = 7c and x[ = Ic, C = Using the relation C = 

we obtain Ac = Sic ^ 2.9 cj^^/^p^/^. We see that the capture amplitude Ac virtually coincides 
with the maximum amphtude of adiabatic oscillations Amax (HOj) . Consequently, it can be 
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assumed that, having passed through the force- free surface and having been deflected to the 
distance A^, the charged particle will be captured and subsequently will begin to execute 
adiabatic oscillations, which we investigated in detail above. For oscillations with the initial 
amplitude Ac, the time Tp reaches its lower limit, Tp ~ (10^^ — 10^')rd. 

For completeness, let us find the frequency of nonrelativistic oscillations u. Without 
providing specific calculations, which are fairly cumbersome, we will write out the final result: 



where rjfs is defined by Eq. (jlj) of the force-free surface. We will also give an expression for 
the radius of curvature of the field lines: 



If p is required to be calculated at some point on the force-free surface, then, naturally, we 
should set r = Tffg in this formula. 

In the entire previous discussion, we generally paid absolutely no attention on precisely 
which particle, an electron or a positron, oscillates near the force-free surface. All of the 
expressions derived above suggested that the oscillating particle had a positive charge. Let us 
turn to Eq. (14T|) for the frequency of nonrelativistic oscillations. In the cases where tu^ > 0, 
positron oscillations with frequency uj take place. If, alternatively, tu^ < 0, then electron 
oscillations take place. Obviously, the frequency of nonrelativistic electron oscillations in this 
case is y—iJ^. The sign of the oscillating particle charge may be said to coincide with that of 

As follows from Eq. fKTl) . the frequency of nonrelativistic oscillations becomes zero only 
if Q = 7T /2 or 6' = 7i/2. The former and the latter equalities correspond to the equator and 
magnetic equator of the neutron star, respectively. The equator and the magnetic equator 
separate the force-free surface into regions in each of which the oscillations of particles only 
with the same sign can take place. When passing through the line formed by the nonisolated 
equilibrium positions, the sign of the particle charge is reversed. The simultaneous fulfillment 
of the equalities 6 = 7c/2 and 6' = 7i/2 defines a straight line at each point of which = as 
well. However, the signs of the oscillating charges coincide on the open sheets of the force-free 
surface that are separated by this straight line. If < 7i"/2, then positrons accumulate in 
the vaulted parts of the force-free surface, while electrons accumulate on its open sheets. The 
signs of the particle charges are shown in Fig. 1. Note that u"^ depends on the wave number k 
and, hence, on its sign. If 6^ > 7r/2, then it is formally necessary to consider the case where 
the angle between the axes of an oblique rotator is vr — 6'^ but, at the same time, to change 
the sign of the angular velocity Q. Consequently, in this case, the signs of the charges will 
be reversed and electrons and positrons will accumulate in the vaulted parts and on the open 
sheets, respectively. 




(41) 



P 



r(l + 3cos2 
3sin^'(l + cos2 9')' 



(42) 



19 



7 The trajectories of particles on the force-free surface 

We saw that the motion of a particle captured by the force-free surface is the sum of the drift 
motion of some guiding center over the force-free surface and the oscillatory motion of the 
particle about this guiding center. Let us find the exact trajectories of the guiding center in 
the rotating frame of reference. 

Assuming that Vq = "^^'ib + v^, it is easy to verify that Vq ■ = 0. We see that 
the velocity vector Vq of the guiding center is always orthogonal to the effective electric field 
E'^-'^-^. Let us find the equation of this surface at each point of which the vector E^-^-^ is directed 
along the normal to the surface in question. If there existed such a function ^ that E^-^-^ = V^, 
then the set of surfaces to which the vector field E^-^-^ is orthogonal would be specified by the 
equipotential surfaces of the function C,- Directly integrating the equation E^-'^-^ = and 
taking into account Eqs. (fT^ . we can make sure that such a function exists: 



The potential ^ is defined to an arbitrary real constant that we do not write for short. The 
equation of the surfaces to which the vector field E*^-^-^ is orthogonal is ^ = €, where <t = const 
is an arbitrary real number. 

Finding the trajectories of the guiding center is reduced to finding the intersection of the 
force-free surface (jlj) and the set of equipotential surfaces ^ = €. Let us introduce a function 



where rffs{0, (p) is defined by Eq. (j3]) of the force-free surface and the potential is defined by 
Eq. The sought-for trajectories are then implicitly specified by the equation C(^)V^) = 

where, as above, €. is an arbitrary real number at which, of course, the equation has solutions. 

Figure 2 shows phase portraits of the trajectories in angular coordinates 9 and (p on 
the surface in the rotating frame of reference for oblique and orthogonal rotators. In what 
follows, we assume that the azimuthal angle corresponding to the magnetic axis is zero. We 
see that there exist no extended trajectories starting and ending at points on the magnetic 
equator. This differs from the result by Jackson [16]. He argued that some of the trajectories 
on the force-free surface of an uncharged orthogonal rotator ended at the magnetic equator; 
as a result, a reverse drift motion to the stellar surface that sweeps out the captured particles 
exists. This could be an obstacle to the formation of a force-free magnetosphere. However, 
Jackson assumed that the particle velocity at some point of the force-free surface was equal to 
the projection of the drift velocity onto the tangent plane at this point. Above, we ascertained 
that this is not the case. 

A 3D image of the trajectories of the guiding center on the force-free surface for an 
oblique rotator is presented in Fig. 3. For convenience, let us introduce the potentials of 




(43) 



R 

km'' 
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isolated equilibrium positions 

^ l(lTcos^^ 
^3(3±cos^„)i/2 

and the surface potential (r = (2/3) cos6'm. In the dome-shaped parts of the force-free surface, 
the trajectories are in the shape of closed loops surrounding the equilibrium position. The 
potential corresponds to the equilibrium position itself, while the interval of potentials 
C_R < C < C- corresponds to the set of loops. The equilibrium positions at the equator and the 
magnetic equator are characterized by the surface potential (r. Two different regions can be 
distinguished on each of the two open sheets in the phase portrait. The first region is similar to 
the region on the vaults and also consists of closed loops surrounding the equilibrium position. 
The potential C+ corresponds to the equilibrium position, while the interval of potentials C+ < 
C < corresponds to the set of loops. However, the loops themselves become increasingly 
elongated as ( approaches zero, tending to go to infinity. The potential C = corresponds 
to the separatrix. The separatrix is an open trajectory that goes to infinity but, at the same 
time, encloses the equilibrium position. The second region on each of the open sheets is formed 
by the set of trajectories with a potential < ( < (r. None of these trajectories lies entirely 
on the same sheet. As we see from Fig. 3, all trajectories pass from one sheet to the other, 
crossing the straight line 6 = 6' = 71/2 and being closed around the neutron star. As C tends 
to (r, these trajectories come increasingly close to the magnetic equator. 

The direction of the trajectories can be easily determined using Eq. (l20l) . In the half- 
space cos6 > 0, the motion along the trajectories in the vaulted part of the force-free surface 
and on the open sheet inside the separatrix = is clockwise when viewed from outside (it is 
implied that the line of sight crossing the force-free surface is directed along the radius toward 
the stellar center). In the half-space cos 6 < 0, the motion along the trajectories in the vaulted 
part of the force-free surface and on the open sheet inside the separatrix is counterclockwise 
when viewed from outside. On the open sheets outside the separatrices (for the trajectories 
crossing the straight line 6 = 9' = 7r/2), the motion is from the half-space simp < to the 
half-space sin (f > when cos 6 > and from the half-space sin > to the half-space sin if < 
when cos^ < 0. All of the aforesaid refers to the case of 6^ < tt/2. If the inclination of the 
rotator x > 7i"/2, then we should consider a rotator with an inclination ft^ = tt — x, direct its 
polar axis from which all polar angles are measured opposite to the vector fl, and change the 
sign of Q in Eq. (1201) . In Fig. 3, this corresponds to reversing the direction of the vector f2. 
Obviously, the direction of all trajectories will be reversed. 

Thus, all trajectories of the guiding center on the force-free surface are closed and lie in 
a finite region, except the separatrix C = 0- As a charged particle moves along its trajectory, 
the electromagnetic field and, hence, the particle oscillation parameters change. However, the 
guiding center moves with drift velocities of ~ 10~^ and the particle oscillation frequency 
exceeds 10 MHz even in the case of ultrarelativistic oscillations. Therefore, the particle travels 
a distance of the order of several millimeters in the oscillation period, which is much smaller 
than R. As a result, the oscillation parameters of the particle as it moves over the force- free 
surface change adiabatically. It is important to note that the sign of cu^ does not change in the 
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case of motion along the trajectory (see the discussion after Eq. fH2|) ). Hence, if the particle 
fell on the force-free surface, then it can no longer leave it, because, first, the trajectory of 
the guiding center is closed and, second, no instability of the particle oscillations along the 
magnetic field develops due to the constancy of the sign of uP' . 

8 Discussion 

The dynamics of the motion of electrons and positrons in the inner vacuum magnetosphere of 
a neutron star can be imagined as follows. A charged particle produced far from the force-free 
surface reaches a relativistic velocity in a time T^ei ~ 10^^^ — 10^^^ s and passes to a quasi- 
stationary regime of motion in a time T^t ~ 10^^ s, traversing a distance of ~ 10 — 100 cm. The 
particle Lorentz factor 70 ~ 10^ — 10* (ITO!) is then completely determined by the condition of 
balance between the power of the forces of an accelerating electric field with E\\/ ~ 10^^—10"^ 
(here, we restore the dimensions of quantities) and the intensity of curvature radiation. In 
this case, the particle moves virtually along the magnetic field line, because the electric drift 
velocity is Ve/c ~ 10^'^ and the centrifugal drift velocity is even lower. During the motion 
along the trajectory, the radius of curvature p and the longitudinal electric field E\\ change 
slowly, leading to an adjustment of the particle Lorentz factor 79. The adjustment time is 
fairly short, tq ~ 10^^'^ — 10^'' s, and the particle travels a distance of ~ 1 cm — 100 m in 
this time; its upper limit is reached near the force-free surface. Since this distance is much 
smaller than the stellar radius, the particle Lorentz factor is, in fact, determined by the particle 
coordinates. As a charged particle approaches the force-free surface, the quasi-stationarity 
condition is violated and the particle passes through the force-free surface with a Lorentz 
factor at the crossing point C^ax ~ 10^ — 10^. Once the particle has crossed the force-free 
surface, it is deflected to a distance A^. ~ A^ax ~ 10 — 100 m fllUl) . Subsequently, the particle 
begins to execute adiabatic ultrarelativistic oscillations with a frequency i/ ~ 10 — 100 MHz. 
The oscillations decay due to the energy losses through radiation and their frequency increases. 
The oscillation energy decreases initially linearly and, in a time Tp ~ 10^® — 10~^ s (see ( 155]) ). 
a power-law decay through curvature losses begins. On time scales of Tcutv ~ ~ 10 — 1000 s, 
the decay becomes exponential with the time constant r^; bremsstrahlung begins to make a 
major contribution to the energy losses. When the decay regime is changed, the particle has 
a Lorentz factor Ccurv ~ 10"^ (l36l) and the oscillation amplitude is A^urv ~ 1 ni. Subsequently, 
the ultrarelativistic particle oscillations continue to decay exponentially and the oscillations 
become nonrelativistic and harmonic with a frequency z/ ~ 1 — 10 GHz when an amplitude 
Inro ~ 1 cm is reached. Simultaneously with the oscillatory motion, the particle executes a 
regular motion over the force-free surface (see Fig. 3); the trajectory of the guiding center in 
the rotating frame of reference is generally closed and its velocity is of the order of the drift 
velocity. 

Let us estimate the accumulation rate of an electron-positron plasma on the force-free 
surface. For this purpose, we considered a situation where the charged particles, electrons 
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and positrons, moved in a given electromagnetic field while acting as test charges. This may 
be assumed only as long as the charge density in the magnetosphere is much lower than the 
Goldreich- Julian density pcj = — B ■ r2/27rc. However, since the pair production rate in the 
magnetosphere is constant, particles increasingly accumulate with time on the force-free surface. 
Consequently, regions will appear near the force-free surface in the vacuum magnetosphere in 
some finite time in which the charge density will now be comparable to pcj- In the long run, 
the entire magnetosphere will be filled with plasma. 

Let us investigate the initial stage of the rearrangement of the magnetosphere from a 
vacuum state to a state filled with plasma. Let us choose an arbitrary point fq of the force-free 
surface and consider some small area on this surface containing the point in question. Since 
the sizes of the area are assumed to be small compared to i?, the area itself is virtually flat. 
Charged particles accumulate with time on this area to form a symmetric charge layer. We will 
consider the plasma accumulation only on fairly short time scales, when the thickness of the 
forming charge layer is much smaller than R. In this case, the plasma density differs noticeably 
from zero only near the force-free surface; therefore, there is as yet no global rearrangement of 
the magnetosphere. For this reason, Eqs. ([T]) and ([2]) specify the external electromagnetic field 
in which the charge layer is located. Let us find the intrinsic electric field of the layer at some 
distance z from the force-free surface: 

z 

El = Ana j pe{z') dz', (44) 



where pe{z') is the volume charge density at distance z' from the force- free surface and the 
distance z does not exceed h — the distance from the force-free surface to the layer boundary. 
Here, it is important that the charge layer is symmetric, i.e., pe{—z') = Pe{z'). To ensure an 
equilibrium state of this layer, it is necessary that each charged particle belonging to the layer 
undergo no acceleration along the magnetic field lines. In other words, the external longitudinal 
electric field must be completely compensated for by the projection of the intrinsic electric field 
of the layer (jHj) onto the magnetic field direction. If we introduce the angle ip between the 
normal to the force-free surface and the magnetic field direction, then 

, B V(E B) 

cos ip 



fi|V(E- B 

For an arbitrary /, the following relation then holds: 

I cos ip 



Ana COS Tp I pe{z ) dz = uj I. (45) 







Differentiating both parts of this equality with respect to / yields pe = u'^/Anacos'^ ip- Thus, 
the charge density pe in the layer does not depend on the distance from the force-free surface 
at all. Recall that we measure the charge density in units of e/X^. It is determined only by the 
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square of the frequency of nonrelativistic oscillations u;^ and the angle tp, whose values are taken 
at point To around which we consider the small area of the force-free surface. The sign of the 
charge density pe coincides with that of u"^, in agreement with our conclusions after Eq. 
In general, pe does not coincide with pcj- This stems from the fact that if the charge density 
in the layer were equal to the Goldreich- Julian density everywhere, then all particles in the 
layer would undergo complete corotation and, hence, would be at rest in the rotating frame of 
reference. In the layer, only the longitudinal electric field is zero but its transverse component 
is nonzero, because the intrinsic electric field of the layer does not necessarily lead to a nulling 
of the transverse component of the external field. For example, since the intrinsic electric field 
of the charge layer at point Tq is zero, the total electric field has only the transverse component, 
which causes the charged particles to drift over the force-free surface along the trajectories 
shown in Fig. 3. Hence, the forming charge layer is composed of differentially flowing currents 
on the force-free surface. 

Thus, in the course of time, a charge layer whose thickness his a function of point fq and 
time t elapsed since the beginning of magnetosphere filling with plasma is formed near the force- 
free surface. It is important to note that only the charges of the same sign as that of the density 
Pe itself contribute to the charge density pe at each point of the force-free surface. Therefore, the 
charges in the layer are completely separated at the initial filling stage in accordance with the 
conclusions after Eq. fH2]) . Consequently, the particle number density is equal in magnitude 
to the plasma density: = \pe\- At each point of the force-free surface, let us introduce 
the particle current density j = dN/dSdt, which coincides with the magnitude of the electric 
current density in dimensionless units. The equation for determining the thickness of the charge 
layer is then dh/dt = j/ng. It can be used to estimate the time it takes for the entire vacuum 
magnetosphere to be filled with an electron-positron plasma. When the thickness of the charge 
layer h becomes comparable to the characteristic size of the inner magnetosphere, i.e., the 
stellar radius R, a significant rearrangement of the entire magnetosphere will take place. The 
filling time Tf is then tj = Rncj/j- Here, we took into account the fact that the plasma 
density in the layer is approximately equal to the Goldreich- Julian density, ne ~ ncj- The flux 
density of the particles falling on the force-free surface is determined by the pair production 
mechanism in the magnetosphere. The gamma-ray photons producing electron-positron pairs 
in the magnetospheric magnetic field can fall into the magnetosphere from outside as a cosmic 
background and can be produced from soft thermal photons emitted by the stellar surface 
by Compton-scattered energetic particles. The main question that arises here is how many 
pairs can be produced by one photon with an energy above 1 MeV. Since the electrons and 
positrons being produced are rapidly accelerated to energies 70 (ITU]) , a chain multiplication of 
pairs takes place. The number of produced pairs per primary photon depends exponentially on 
the characteristic size of the inner magnetosphere, exp/x, where p = R/i. Here, i is the total 
mean free path of an energetic particle with respect to the generation of a curvature photon 
and its mean free path with respect to the production of a pair. The mean free path i is known 
from an investigation of the stationary plasma generation processes in a polar magnetosphere 
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17l | to be ~ 100 m. This means that the muhiphcation factor can reach enormous values, of 
the order of exp 100 ^ 10^^. However, since there is an exponential factor, this question in 
the case of a vacuum magnetosphere requires a separate careful analysis, which is beyond the 
scope of this paper. A large pair multiplication factor can compensate for the smallness of 
the cosmic background photon flux. According to observations Q, Q] , the diffuse background 
jph of Galactic photons with energies above 1 MeV is ~ 10^'^ cm^^s^^. Therefore, the ratio 
jph/cncj ~ 10~^^ is very small, but it does not prevent the pulsar magnetosphere from being 
rapidly filled with an electron-positron plasma at a sufficiently large pair multiplication factor. 
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Figure 1: Sections of the force- free surface by the planes ip — ipm = {0, vr} (left) and ip — ipm = 
{7r/2, 37r/2} (right) for inclinations (from the top downward) of 0, 7r/6, 7r/3, and 7r/2. The signs 
of the accumulating charges are marked. The sections pcj = are also shown for comparison. 
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Figure 2: Phase portraits of the charged particle trajectories on the force-free surface in 6 and ip 
coordinates in the rotating frame of reference for an obhque rotator with an inchnation of 7r/3 
(left) and an orthogonal rotator (right). 
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Figure 3: Charged particle trajectories on the force-free surface in the rotating frame of reference 
for an obhque rotator with an inchnation of 7r/3. The heavy hnes indicate the equator, the 
magnetic equator, and the separatrices. 
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